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Differential phase extraction in an atom gradiometer 
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We present here a method for the extraction of the differential phase of an atom 
gradiometer that exploits the correlation of the vibration signal measured by an 
auxiliary classical sensor, such as a seismometer or an accelerometer. We show 
that sensitivities close to the quantum projection noise limit can be reached, even 
when the vibration noise induces phase fluctuations larger than 2tt. This method 
doesn’t require the correlation between the atomic and classical signals to be perfect 
and allows for an exact determination of the differential phase, with no bias. It 
can also be applied to other configurations of differential interferometers, such as 
for instance gyrometers, conjugate interferometers for the measurement of the fine 
structure constant, or differential accelerometers for tests of the equivalence principle 
or detection of gravitational waves. 


I. INTRODUCTION 

Atom interferometers have demonstrated performances comparable or better than state- 
of-the-art classical instruments, both in terms of sensitivity and accuracy, and hnd applica¬ 
tions in various helds, from fundamental physics to geophysics and navigation. In particular, 
differential atom interferometers allow for an improved determination of the quantity to be 
measured, as their mode of operation rejects common mode noise sources and systematic e 
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fects. They are used for the measurement of gravity gradients [l, 2| and rotation rates M. 




7| and the ratio h/m [8|, l9|, and for tests of the universality of 
10l-ll2| and are expected to hnd applications in the detection of gravitational waves 
16|. In most cases, the huctuations of the common phase, in general dominated by the 


effect of parasitic vibrations, wash out completely the visibility of single interferometers, but 
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the correlation between the two ontput signals can be used to recover the information on 
the difference between the two interferometer phases. Various techniques have been demon¬ 
strated to extract the differential phase. A method based on ellipse fits was first used in jl^ . 
but was shown to introduce signihcant biases, that are errors in the determination of this 


differential phase. Methods based on Bayesian estimators are more accurate 


18 


20|, but 


require an a priori model of the interferometer phase noise. In a recent study, it was shown 
that the use of a simultaneous third measurement allows for the retrieval of the differential 
phase(s) with negligible bias(es) using a robust three dimensional ht 21|. 

In this paper, we present an alternative method, that uses the correlation of the inter¬ 
ferometer phases with an estimate of the vibration phase provided by an auxiliary classical 
sensor. It allows to recover the differential phase of interferometers operated in a gradiome- 
ter conhguration, without any biasing or an a priori knowledge of the noise distribution. 
We show in particular that sensitivities close the quantum projection noise limit can be 
obtained. 


II. PRINCIPLE OF THE METHOD 


The output signals of the interferometers are given by the transition probabilities Pi and 

P2- 

Cl 

Pi = + ~cos(0m + 01 + 4>ll) 

c 

P 2 = A 2 + —cos(0m + 02 + 02 z) 


where 0m is the phase shift induced by the vibrations (of the mirror used to retroreflect the 
lasers of the interferometer), 0j are the two interferometer phase shifts (with contributions 
from inertial effects, such as gravity acceleration, and non-inertial effects, such as lights 
shifts) and (pu are controlled phases shifts applied onto the phase of the interferometers. 
0j/ can be common to the two interferometers (such as when applying to the laser phase 
difference either a frequency chirp for Doppler shift compensation or a phase jump), or 
not. In particular, such differential phase shifts can easily be set in the conhguration of a 
gradiometer, taking advantage of the spatial separation between the interferometers. This 
was realized for instance by applying a magnetic pulse at one of the two clouds [i|,Q,orby 
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inducing a common frequency jump onto the interferometer lasers 23|, the latter technique 
providing in principle a better control of the applied differential phase. 

Without loss of generallity, we will take = 0. This can easily be realized, for 

instance by applying a phase jump before the last pulse of the interferometer. This leads to 


P, = A 


Pi = A^ 

C2 


+ ^COs(0m) 
COs{(j)m + (j)d) 


where = 02 ~ 0i + 02 Z ~ ‘Pii is the differential phase between the two interferometers. 

An estimate 0^ of the phase shift induced by the vibration of the mirror 0^ can be 
determined thanks to the measurement of an auxiliary sensor, such as a seismometer 2^ or 


a mechanical accelerometer 


25 


26j, or with an additional interferometer as 


in 1^]. 


When 


using a classical sensor, a rigid link with the mirror is required for an optimal correlation 
between 0^ and 0s, and some signal processing can be used to flatten the transfer function 


of the sensor 


24|. 


This correlation can be exploited to determine the phase of the interferometer, even 
in the presence of large vibration noise that washes out the fringes at the output of the 
interferometer. One can exploit for instance the fringe fitting technique described in 2^ , that 


is based on the recovery of interferometer fringes when plotting the transition probability 
versus 0s. Vibration noise provides there a random sampling of the interferometer phase 
that allows the fringe pattern to be scanned. 

We show here that this technique can be extended to the case of a differential inter¬ 
ferometer, by reconstructing the fringes for both interferometers independently, and htting 
them as a function of 0s. Though the dispersion of the phases of the two interferometers, as 
determined by the fits, is linked to the quality of the correlation between real and measured 
mirror vibrations (and might thus exhibit relatively large fluctuations when this correla¬ 
tion is poor), the differential phase can potentially be much better determined, because the 
fluctuations of the fitted phases of the two interferometers are also correlated. 
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III. SIMULATION 

To simulate the signals at the output of the interferometer, we start by randomly drawing 
the values of (pm in a gaussian distribution with a standard deviation of In order 

to account for an imperfect correlation between the mirror motion and the sensor signal, 
we randomly draw the difference between the corresponding phases 5p = pm — Ps, with 
a standard deviation of ascj,- The two transition probabilities are then calculated, with 
Ai = A 2 = 0.5 and Ui = U 2 = 1. We then account for the influence of detection noise, that 
we take as quantum projection noise limited with a number of detected atoms of Nat = 10® 
for each interferometer. We thus add to Pi noise contributions (5p., randomly drawn in 
gaussian distributions with standard deviations ap^ = Ppl — Pp/Nat- 


IV. RESULTS 


First, we generate 10® pairs (Pi, P 2 ) with = 3 rad, = 0.3 rad and pd = 0. Such a 
level of vibration noise and degree of correlation between mirror vibrations and sensor signal 


corresponds to the case studied in 


28l |. where an atom gravimeter operating in a urban 


environment was directly put on the ground and the classical sensor was a seismometer. We 
then group the data by sets of 100 trials and perform 10^ consecutive hts of P* versus ps-, from 
which extract 10^ values of the htted phases pij, and of their difference pdj = p 2 j — Pij- 
Figure [H displays the Allan standard deviation of the difference of the htted phases j- 
(open circles), compared with the htted phase of one of the interferometers (open squares) 
and with the vibration phase noise (full squares). While the correlation with the data of the 
classical sensor allows to gain a factor of about 8 on the sensitivity of a single interferometer, 
the gain on the diherential phase noise is much larger, of about 1700. The thin line represents 
the quantum projection noise limit of aqpN = V^/VNat, and only lies a factor 1.25 below. 

We then investigate the efficiency of this rejection as a function of the value of the 
diherential phase, by plotting the ratio ^ /aqppi as a function of the diherential phase 

pd, for a hxed value of = 3 rad, and diherent amplitudes of asif,- The results are 
displayed on hgure [2J We observe an optimal phase sensitivity when pd = 0[7r], close to 
the quantum projection noise limit. We observe a rapid degradation as the value of pd 
starts deviating from these values, that increases with as,p- In particular, the correlation 
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FIG. 1: Allan standard deviation of the differential phase (open circles), compared to Allan stan¬ 
dard deviations of the phase of a single interferometer (open squares) and of the vibration noise 
(black squares). The line represents the quantum projection noise limit. Parameters are = 3 
rad, (T 50 = 0.3 rad and = 0. 

between the two fitted phases is lost when operating the two interferometers in qnadratnre. 
There, the sensitivity is minimal and comparable to the sensitivity one wonld obtain with 
independent nncorrelated interferometers. 

In a second series of simnlations, we vary the amplitnde of as^p, keeping 0^ = 0 and 
^<t>m = 3 rad. The results displayed on hgure E] show that the sensitivity in the differential 
phase remains close to the QPN limit, within a factor of 2, over a large range of ascj, of 1 
radian. This shows that the retrieval of the differential phase is still possible, even when 
the correlation between the classical signal and the mirror motion is far from being perfect 
and the level of nncorrelated noise is high. Degradation of the correlation can arise from 
non linearity of the sensor in the frequency range where the vibration noise is important 
(from 0.1 to 100 Hz), or from the intrinsic noise of the classical sensor. In the example 
above, a level of intrinsic noise of the classical sensor lower than the level of vibration noise 
by a factor of only 3 would still allow reaching a sensitivity as good as twice the quantum 















6 



<t>d (rad) 


b) 



0.00 0.02 
<l>d (rad) 



FIG. 2: Sensitivity in the differential phase extraction, normalized to the quantum projection noise 
limit as a function of the value of the differential phase. A zoom on the curve for small values of 
the differential phase, below 50 mrad, is displayed on the right. 

projection noise limit. This shows that, depending on the parameters, mid class sensors, 
which are more compact and less expensive, can in practice still be highly benehciah 

We then repeat these simulations for increasing values of 0^, ranging from 10 mrad to 
200 mrad. Figure 0] shows the results, that illustrate that the trend in the degradation of 
the sensitivity versus 6(j) increases with increasing values of 4>d, which was already observed 
in hgure [H 

Finally, we vary the amplitude of keeping the quality of the correlation the same, 
which means that we keep the ratio constant, equal to 10. Figure [5] shows an 

increase by less than a factor 2 with below an amplitude of 10 rad. Note that for such 
a value, = 1 rad which corresponds to a situation where the fringes obtained when 
plotting Pi versus (prn are hardly visible. Even in that case, the sensitivity is only twice the 
QPN limit, and the rejection efficiency of vibration noise is as large as a factor 4000. For 
larger vibration noise, the degradation of the sensitivity gets significant. Also, the number 
of points to be htted needs to be increased for the ht to converge to the correct value. We 
obtain stable fits up = 20 rad with a number of points to fit of 1000. 

Of importance, we checked that, whatever the noise parameters we choose, the extraction 
of the differential phase is free from any biasing: the mean value of the fitted differential phase 
is always found equal to the differential phase. This will hold in a practical implementation 

















7 



FIG. 3: Sensitivity in the differential phase extraction, normalized to the quantum projection noise 
limit, as a function of the standard deviation of 64>, the difference between real and measured mirror 
vibrations. The differential phase is zero. 


in particular when the mean value of and cj)m are equal (and for instance equal to zero 
in the absence of a lon g te rm motion of the mirror). This equality is verihed if using a low 


noise seismometer 


24 


, whose output signal is a velocity signal, but not necessarily if 


using a DC accelerometer, such as in j2^. In the latter case, a non zero DC value leads to 
an offset in the htted phase, that can exhibit long term drifts. Nevertheless, this offset being 
common to both interferometers, it will cancel out when calculating the differential phase. 
Finally, we checked that the results presented here do not depend on the values of the offsets 
Ai and A 2 on the transition probabilities, nor on the values of Ci and C 2 (if accounting for 
the corresponding degradation of the QPN limited sensitivity). 


V. CONCLUSION 

We presented here a technique for extracting the differential phase in an atom gradiometer 
in the presence of large vibration noise, that washes out completely the visibility of the 
fringes of the individual interferometers. It relies on the exploitation of the auxiliary signal 
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FIG. 4: Sensitivity in the differential phase extraction, normalized to the quantum projection noise 
limit, as a function of the standard deviation of 6(j), the difference between real and measured mirror 
vibrations. The differential phase is varied between 10 mrad and 200 mrad. 


provided by a classical sensor, that measures the motion of the mirror that retroreflects the 
interferometer lasers. We show that sensitivities close to the QPN limit can be reached, 
provided that the value of the differential phase modulo vr is zero. This corresponds to 
in-phase or in counter-phase operation of the two interferometers, which can be set either 
by changing the scale factor of the interferometers, or by applying a perturbation that is 
not common mode to the two interferometers. The technique can also be used for other 
differential interferometers, such as gyrometers based on counterpropagating atomic sources 
4, Isl], or conjugate interferometers for the determination of the ratio h/m 29|. 

This technique is also of interest for differential accelerometers that use two different 
atomic species in order to perform tests of the equivalence principle. In that case though, 
the accelerometer scale factors may differ due to the difference in the effective wave-vectors 
of the interferometer lasers. Vibration phases, even if proportional, are thus different for 
each interferometer 


20| . As an alternative, the scale factor can be matched 


interferometer pulse sequence slightly different for the two accelerometers 


Dy making the 


19 1. Though the 
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FIG. 5: Sensitivity in the differential phase extraction, normalized to the quantum projection 
noise limit, as a function of the mirror vibration noise The ratio between and a 5 ^ is kept 
constant, equal to 10. The differential phase is null. 


response of the two interferometers to DC acceleration are then identical, the difference 
in the pulse sequence makes the transfer function not identical for the two accelerometers 
and thus their response to vibration noise is different 
between the htted phases is then reduced, a recent study 
technique compares favorably with other methods, and allows for a more robust extraction 
of the differential phase, free of any bias, and with a potentially excellent sensitivity. 


30i. Even though the correlation 
301 1 demonstrates that the present 
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